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Abstract 



Starting from the QCD Lagrangian we derive the effective action for 
massive quark and antiquark at large distances, corresponding to the minimal 
area low of the Wilson loop. The path integral method is used to quantize the 
system and the spectrum is obtained with asymptotically linear Regge trajec- 
tories. Two dynamical regimes distinguished by the string energy-momentum 
distribution are found: at large orbital excitations (/ ^ 1) the system behaves 
as a string and yields the Regge slope of while at small / one obtains a 
potential-like regime for relativistic or nonrelativistic system . The problem of 
relative time is clarified. It is shown that in the valence quark approximation 
one can reduce the initial four-dimensional dynamics to the three-dimensional 
one. 

The limiting case of a pure string (without quark kinetic terms) is studied 
and the spectrum of the straight-line string is obtained. 
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1 Introduction 



This is the first paper of the presumed series devoted to the quantum dynamics 
of the quark- antiquark system at large distances. Our starting point is the 
formalism of vacuum correlators, developed previously [1] (for a review see [2]). 
It allows one to represent the gauge-invariant Green's function of the quark- 
antiquark system in a form, where all dynamics is contained in the averaged 
Wilson loop operator. 

We simplify our problem by disregarding effects due to the quark spins and 
additional quark loops (sea quarks) having in mind to come back to it in later 
papers. 

In this way our simplified problem is that of a scalar quark and antiquark 
without additional quark loops in the confining background field. 

One of the most important point of the paper is to identify an explicit mech- 
anism creating the QCD string and to find the properties of the latter. 

Even the question: what is it, the QCD string? is not trivial. Usually it 
is associated with the Nambu-Goto string action (for review see [3]) describing 
the open string with active degrees of freedom all along the string including 
the ends, where massless quarks are presumed to be described by the proper 
boundary conditions. This standard picture is plagued by unphysical features 
in 4 dimensions and needs o? = 26 or supersymmetric extentions to be consistent 
[3]. 

The picture which emerges in our paper and based on the QCD Lagrangian 
and the vacuum correlator method essentially differs from the standard one. 
First, the string appearing in the qq system has a world-sheet coinciding with 
the surface appearing in the area law of the Wilson loop. This area law is a 
natural consequence of the cluster expansion for the Wilson loop [4], and the 
surface bounded by the Wilson contour appears also naturally in the formalism. 

An important point is what kind of surface appears in the area law of the 
Wilson loop ? Since the whole cumulant series does not depend on the shape 
of the surface and has no degrees of freedom on it , we come to the conclusion 
that this should be the minimal surface which enters the area law and therefore 
defines also the shape of the string, connecting quark and antiquark. We stress 
however that this is an assumption we make here |^, and it is desirable to obtain 

^This assumption is equivalent to the uniform convergence of the cluster series for the Wilson loop [4], since 
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the behaviour of the Wilson loop in QCD directly. 

The QCD area law appears only asymptotically at large distances, R ^ T^, 
where Tg is the so-called gluon correlation length, which enters as a scale for the 
vacuum correlators [5,1,2]. Recent Monte-Carlo calculations give the estimate 
Tg « 0.2 0.25/m [6] and therefore we can apply our formalism even to the 
ground state mesons (see [2] for a discussion). We approximate the minimal 
area law surface as the world sheet surface of the straight line string connecting 
proper positions of quark and antiquark. We shall discuss in the next sections 
the validity of the straight-line approximation. 

We consider the string tension as being renormalized and disregard all cor- 
rections including perturbative gluon exchanges and the creation of additional 
quark loops, since our aim here is to concentrate on the main dynamical in- 
gredient - the interaction between quarks corresponding to the minimal area 
law. 

Therefore in this approximation the minimal string may rotate (orbital ex- 
citations), and oscillate longitudinally (streching and expanding). The latter 
type of motion is not possible for the standard Nambu-Goto string. 

After these introductory words about the definition of the QCD string, we 
can outline the purpose and the plan of the paper. 

We shall quantize the system consisting of a massive or massless quark and 
antiquark connected by the minimal string. 

We shall extensively use the method exploited recently by one of us [2,7] 
for this purpose. The method contains five steps. (1) First, one represents the 
quark- antiquark Green's function using the Feynman-Schwinger form [1] as a 
path integral over trajectories of q and q parametrized by their proper time 
parameters. (2) Second, all gluonic nonperturbative field contained in Wilson 
loops, is replaces by the minimal area law. (3) Third, the minimal area is 
assumed to be spanned by straight lines connecting a point on the q trajectory 
with another point on the q trajectory, these two points are chosen at the same 
proper time parameter. (4) The proper-time Hamiltonian is introduced to get 
rid of path integration and obtain instead differential equation. (5) Fifth, an 
approximation is made in the method [2,7] to expand the Nambu-Goto string 
term in the action in powers of some numerical parameter (which yields also 
expansion in powers of relative velocity r^). 

Due to this expansion dynamics in relative time is free and one integrates 
it out, regaining in the end three dimensional (relativistic) dynamics with dy- 

any term in the series provides the area law for large enough loop. 
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namically generated constituent mass. 

This approximation [1,2,7] appeared to be simple and practical, yielding 
masses and Regge trajectories in terms of only one parameter - the string 
tension. In this way light mesons, heavy-light mesons [8], baryons [9] and 
glueballs [10] have been considered. The accuracy of the approximation for 
masses was estimated as ~ 10%. 

As a result [7] the slope of Regge trajectories was obtained to be (8cr)"^ the 
same as in the potential models [11] in contrast to the Nambu-Goto string slope 
of {27Ta)-\ 

In the present paper we abandon this approximation and treat the effective 
action of the system exactly, using the formalism of auxiliary fields to get rid of 
the square root term [12,3]. As a result we shall not only improve the accuracy 
of the approximation but shall obtain a qualitatively different dynamics, which 
has not been present before in the approximation [2,7]. 

We shall find that there are two regimes of qq dynamics distingushed by the 
energy- momentum distribution of gluon fields, — for large orbital momenta 
/ the resulting spectrum in the leading approximation coincides with that of 
a pure string with a slope l/27rcr , while for low values of / (depending on 
the masses of quarks) the dynamics is described by the relativistic potential- 
like approach, close to the one obtained previously [8,2]. For the heavy quark 
system the potential picture is valid in a large region of / and it joins smoothly 
the string picture for very large /. 

The limiting case of pure strings (without kinetic terms of quarks) is con- 
sidered in detail Q. The straight-line Regge-trajectories, corresponding to the 
spectrum of the rotating straight-line string is obtained. This is in an agree- 
ment with the results obtained in [14]. Note that in our approach the string 
picture has not been assumed but it was derived from the QCD-Lagrangian 
under rather general assumptions. 

The important issue of relativistic dynamical systems is the role of the rela- 
tive time of q and q. To make our straight-line approximation for the minimal 
surface selfconsistent we have to integrate over the class of paths without time 
backtracking of quarks (in the c.m.s.). This constraint corresponds to the va- 
lence quark approximation to the problem and allows one to rigorously reduce 
the initial four-dimensional (4D) dynamics to the three-dimensional (3D) one. 
As a result the quadratic 4D kinetic terms are transformed into 3D square-root- 
type terms. We have the Lorentz invariant effective action with the constraint, 

short version of these resuhs is reported in ref. [13] 
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which can be formulated in a Lorentz covariant form. So the calculation of the 
spectrum can be performed in an arbitrary frame and for the sake of convenience 
we work in the meson rest system. 

The plan of the paper is the following. The Feynman-Schwinger represen- 
tation (FSR) for the qq - Green's function is given and the effective action for 
qq-system is obtained in Sec. 2. In section 3 we obtain the approximate explicit 
expression for Wilson loop at large distances. In section 4 in the valence quarks 
approximation we integrate out time-components, so that we reduce the 4D dy- 
namics to the three-dimensional one. Gaussian representation for the obtained 
action is formulated in Sec. 5. The method of auxiliary fields [3], [12] is used 
in this section to get rid of the square root term, which determines the string 
action. 

The case of a pure straight-line string without quarks at the ends is discussed 
in Section 6. In Section 7 Hamiltonian for the general case of the straight-line 
string with quarks at the ends is obtained and in two limiting cases I « 1 and 
/ ^ 1 the analytic form of the spectrum is established. 

In conclusions we summarize our main results and make comparison to those 
in literature. 

Appendices A,B,C contain technical details of the derivation of formulas in 
the text. 

2 Green's function of the quark-antiquark system inter- 
acting with gluon field. 

In this Section we use the Feynman-Schwinger representation for the quark- 
antiquark Green's function to obtain the effective action for qq system in terms 
of the Wilson loop [1]. 

We start with the initial and final qq states defined on space-like surfaces in 
a gauge-invariant way 



where Fj^, Tout contain a parallel transporter $ and some vertex with definite 
Lorentz structure 7^^, jout- 



'^in{y,y) =u{y)Tin{y,y)u{y) 



(1) 



(2) 



Tin = y)-fin, y) = Pexp{ig A^dz^) 



(3) 
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and the same for Tout] lin may be e.g. 75, 7^ etc. The operator P ensures the 
the ordering of A^{z) along the path Zfj^{t). The Green's function G{xx/yy) 
is obtained by averaging the product ^^^^^^^ over all quark and gluonic fields 
inside the path integral with the usual QCD action 

G{xx/yy) =< x) >^,^ (4) 

The quark degrees of freedom can be easily integrated out with the result 
(for a nonzero flavour channel) 

G[xx/yy) =< trr{x, x)S{x, y)T{y, y)S[y, x) >a (5) 

where S{x^ y) is the quark propagator in the gluonic field A^, and the averaging 
over gluonic fields now includes also the quark determinant (actually a product 
of determinants over all flavours): 

< B >A= J DAe-^^^^B{A)l[det{mi + D{A)) (6) 

i 

For the quark propagator in the gluonic field one can use the Feynman- 
Schwinger representation [1,2,15] 

S{x, y) =< X \ {m + D{A))~'^ \ y >= 
=<x\{m- D{A)) r (ise-^(™'-^'(^)) | y >= (7) 



^ poo 

= (m - D{A)\\^ dsDze-'^^Y.ix.y) 
where the following notations are used 

K = m'.^\llilm. ^. = ^ (8) 

*s(ir, y) = (PEexpig /J J:,,„F^Mt))dt))Hx, y) (9) 

and S^^j, = |;(7/x7j/— 7!^7/i)- The operators ensure the proper ordered insertion 
of operators T^^j^yF^j^y along this path (for a discussion see [15]). 
Insertion of (7) into (5) immediatly gives [8,15] 

G{xx I yy) = 

I dse-''Dzj^ dse-''Dz<^i,^^[x,x){m-D[A)),^^[x,y)x (10) 

x-iin^{y,y){m - b{A))y^Y^{y,x) >a 
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Analysis of the factors (m — D) made in [9] shows that it can be replaced 
by (m + and taken out of the angular brackets, so that we are left with 
the only term defining all the dynamics of the system, both perturbative and 
nonperturbative, which we call W^.- 

The presence of spin degrees of freedom in (7) makes dynamics rather com- 
plicated. For this reason we shall omit the spin dependence and concentrate 
on the confining dynamics without spins. We also neglect the additional quark 
loops, which corresponds to the substitution of Ildet in eq. (6) by 1. This 
is in line with usual quenched approximation and can be justified in the 1/Nc 
expansion. 

The resulting Green's function in this case has the form 

G{xx\y,y) = ds dse'^'^ DzDz <W{C) >a (H) 
where W{C) is the usual Wilson loop operator 

W^{C) = tr^{x,x)^Y,{x,y)^{y,y)^^{y,x) = trPexp[ig j^A^dz^] (12) 

The closed contour C consists of initial and final pieces [x, x], [y, y] and paths 
2;(t), z{t) of the quark and antiquark. 

It is convenient introduce as in [2] the following " center of mass" coordinate 
i?^ and the "relative" coordinate 

i?^ = -^S-jit) + -At)) 

s -\- s s s 
= z{t) - z{t) 

with the boundary condition imposed in the Lorentz covariant way 

i?/z(s) - i?^(0) = Tw^ wX = l 

and r^{s) = 'r^i(O) are also fixed. 

In what follows we shall be interested in the limit T ^ oo to consider asymp- 
totical states of the system. In this hmit our boundary condition corresponds 
to the total momentum of the meson ~ u^. 

3 Evaluation of the Wilson loop at large distances 

In this section we shall obtain approximate expression for < W{C) >a at large 
distances. 
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Using the cluster expansion we can represent < W{C) > ) as [4] 

< WiC) >= exp f: ^ / dail)...daik) < F(1)...F(A;) > (13) 



k=l 



where cumulants are irreducible averages [16,17] and we omit Lorenz indices in 
d(j^^{u{i)) and F^^{u{i)). 

As shown in [2,18], there are three possible regimes for < W{C) > de- 
pending on the relation between the sizes of the loop C and the correlation 
length Tg , defining the decay of cumulants <C F{l)...F{i)...F[j)...F{k) 

If we represent the loop C as a rectangular of time length T and space 
width i?, then we have (omitting perimeter-type terms always coming from the 
quark-mass renormalization and exchanges) [2] ' 

i) < W{C) >^ exp{-(TS),R :>Tg,T:> Tg (14) 

a) < WiC) exp{-T{c2R^ + C4i?^ + ...)), R<.Tg,T:$>Tg (15) 

in) < W{C) >^ ^^^Pi- ^'^' 24n!^^ ^ + ^ ^ ^ ^^^^ 

Here S is the area of the surface inside the contour C, which we take as the 
minimal surface, since < W{C) > does not depend on the shape of the surface 
[17]. 

Thus we see that the area law (14) and the ensuing string dynamics appears 
only as the asymptotic regime for large contours, while for small or narrow 
loops the area law (and the QCD string) is absent and is replaced by a much 
weaker interaction. 

In what follows we concentrate on the regime (14), i.e. we shall consider 
only large loops, R,T ^ Tg. 

We note that in the Monte Carlo calculations this law has been seen in the 
wide region of R,0.1 ^ R ^ Ifm [19] also in the presence of dynamical quarks 
(the determinantal term in (6)) [20]. From those results one could conclude 
that Tg is small enough, so that the regime (14) is dominant for light quark 
systems. 

Independent Monte-Carlo calculations [6] have confirmed that Tg = 0.2 
0.3/m, which should be compared with characteristic quark relative distance 
and orbiting time R Tq ^ Ifm for light quark system [2]. Thus we can 
conclude that the area law (14) is a reasonable approximation for light quark 
systems (even for ground states) and for excited states of heavy quarkonia. 
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The "minimal area law" implies the neglect of gluonic excitations above the 
QCD vacuum, which could lead to an effective integration over the surfaces, 
bounded by the contour C . 

Next step is to construct explicitly the minimal area S in terms of a given 
contour C, defined by quark and antiquark paths z^{t) and z^^it). 

Any surface can be parametrized by the Nambu-Goto form 

^ = lo lo - (^/^%)']'^' (17) 

where w;^(r, (3) are the coordinates of the string world surface, and 



We use the approximation that the minimal surface for given paths z^{t), 
Zij{t) is determined by eq.(17) with w;^ given by straight lines, connecting points 
Zi^{t) and Zi^{t) with the same r, i.e. 

Wi,{T,^) = z^{r)-(5 + z^{r){l-(5) , ^ /? ^ 1 (18) 

where r is defined for both trajectories as 

t i 

T = -T = -T (19 

s s 

and T enters boundary condition. We have no proof that eq. (18) indeed fits 
the minimal area S in all cases; we note that this approximation is valid in 
two limiting cases which are of special interest below: in the case / = one 
can exploit the fiat dynamics of quarks and in the limit / — ^ oo corresponding 
to the dynamics of the string with typical trajectories of the double helycoid 
type, for which the minimal area indeed is formed by the straight-lines. In what 
follows we shall use eqs. (14,17,18) to express < W{C) > in terms of quark 
coordinates. 

We introduce a compact notation 

^ = {r, f3} , gabiO = daw^dbw^, a,b = T,f3 (20) 

and obtain 

S = J Si^d^g. (21) 
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4 The reduction of four-dimensional 

dynamics to the three-dimensional one 

Before the insertion of the concrete form of < W{C) > let us consider the 
integration over zo{t),zo{t) and its physical interpretation (in the meson rest 
system). These integrations (additional with respect to the nonrelativistic case) 
are related to the following phenomena. First, the relativistic quantum theory 
in general is the theory with unconserved number of particles. Even in the case 
of a free particle the backtracking of the time-component which corresponds to 
the creation of additional pairs at the intermediate stage leads to the proper 
relativistic structure of Green's function. Second the fact, that the interaction 
time is smeared out, generates a need for the integration over the relative time 
of particles in the process even in the valence quark approximation. 

We suggest the following approach to the problem. At first step we separate 
out from the initial sum over all paths the class of the trajectories without 
backtracking of time-components, satisfying the condition (in the rest system) 

^^0(71) ^ Q ^^^0(72) ^ Q 
or rewriting it in the explicitlly Lorentz-invariant form 

(P^ dz,{-fi)/d-f,) > , {P^dz^{j2)/dj2) > (23) 

where is the total 4-momentum of the meson. 

This constraint has to be imposed, because for the trajectories with such 
backtracking the minimal surface which appears in < W{C) > obviously can 
not be properly approximated by our straight line anzatz. This procedure 
corresponds to the usual valence quarks approximation and in the language of 
the time ordered diagrams the trajectories with backtracking corresponds to 
the production of qq - mesonic states. 

In Appendix A it is shown, that for the trajectories without time-back- 
tracking the variables ^0(71) , ^0(72) can be transformed into the new nondy- 
namical ones I2i{r) , [J'2{t) without derivative terms by using the parametriza- 
tion for which 

z,, = {t ,z) ,z^ = {t , z) (24) 

The integration over these new variables for the spinless quarks is performed 
effectively by steepest decent method and leads only to the modification of the 
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(22) 



initial kinetic terms K (8) in the quark Green's function in an external field 
(7). The final action becomes independent of the time-components of z^, z^. 

G[xx/yy) = J Dfn D112 Dz Dz exp[-K' - K'] < W{C) >a (25) 

where 

K' + K' = ]^[i-^ + ^,^ir){l+^^ir)}) + i-^^ (26) 
i 2 ^i(rj fi2[T) 

We have introduced the new proper time parameter r , < r < T , and 

T T ^ 

^i(^) = ^^o(^) ' ^2(r) = ^^o(r) (27) 

with the dot standing for the derivative over r. It should be stressed that the 
substitution of eq.(24) into < W{C) > is implied. This condition is equivalent 
to 

i?o(r)=r, ro(r) = (28) 

which corresponds to the usual instant plane Hamiltonian dynamics in the cm. 
system. We also observe in (28) that r plays the role of the time for the meson 
- a common time for both quarks. We shall prove, that in the case of the pure 
straight-line string the condition 

(Pr) = 

directly follows from the dynamics of the string. This result supports the va- 
lidity of our approximation of neglecting (in the straight-line approximation for 
the minimal surface) of backward in time trajectories in the rest system. It is 
important to stress, that this assumption can be not valid for the calculation 
of quantities other than the mass spectrum and for interactions different from 
the confining string-like interaction of quarks, which we have formulated. 

Therefore we have shown that the starting 4D dynamics can be reduced to 
the 3D one. And it should be emphasized that in general this is a nonlocal 3- 
dimensional dynamics due to the double independent integration in W{C) over 
z{t) and z{t), which bindes together z{r) and z{t) with different arguments ri 
and r2. 

In addition we have two additional integrations over Dfii, and Dfi2 which 
replace integrations over time components zq^zq, and one common evolution 
parameter in the action which can be identified with the cm. time. The 
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physical meaning of fii, fi2 can be most easily clarified imposing an external 
e.m. field, e.g. in calculating the magnetic moment. Then it can be shown, 
that fii enter magnetic moment as the constituent mass of a quark. The same 
happens with spin-dependent forces (to be published). Therefore we shall call 
Hi the dynamical mass of the quark i. In the approximation when fi does not 
depend on r, this dynamical mass has been introduced in [2,7]. As we shall see 
below this approximation works reasonably well - accuracy in the determination 
of mass is around 5% ( see e.g. Table 4 of [2]). 



5 Gaussian representation for the effective action of 
quarks and the string. 

Combining the results of the previous section we obtain the total effective ac- 
tion. 

A = K' + K' + ao dr d(]^fd^ (29) 

This action is similar to the one considered in the papers [15] , where the 
straight-line string without transverse excitations has been considered. However 
in our case also kinetic terms of quarks K' + K' are present and the condition 
that the ends off the string move with the velocity of light is not possible. 

A direct procedure of quantization of (29) is difficult due to the square root 
term and we use the auxiliary fields approach [3] to get rid of it. In the Appendix 
B we give a detailed derivation of this procedure , while here we present the 
final gaussian representation of this action 



A = r dr f dP{\(^ + ^) + \,,4T)te + + (30) 



where 



+ ^[w^^ + (o-i>)V^ - 277(i(;r) + T/V^]} 



/i+(r) = ^i(r) + /i2(r) , ^(r) = 



/^i(r) +//2(r) 

and the condition (24) for 2;^(r),z^(r) is assumed. Since only extremal values 
of jJLiij) and //2(t) contributes one can intgrate in our symmetric case over 
restricted class of functions //i(r) = //2(t) = l^ij). Here y{T^j3) > and r]{T,/3) 
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are two auxiliary fields, which should be integrated out in the full path integral 
representation for G: 

G = j DRDrDDDrjD/ie-^ (31) 

In order to prove the equivalence of the dynamics governed by the ac- 
tions (29) and (30) it is enough to show that after proper integration over 
z>(r, /?), 77(r, /?) one returns from eq. (30) back to the initial one (29). 

Actually after the gaussian integration over ry(r, /?) we obtain instead of eq. 
(30) the following expression for G 

G = I DRDrDiiexp[-K' - K']Dux (32) 

where we restore the notation 

For the integral over z>, as it is shown in Appendix A, the explicit measure 
can be constructed and the effective action is determined only by the extremum 
value of i>. 

o\^vl = [wV - [wwf) (33) 
Inserting this expression for 9 into the eq.(32) one recovers our starting action 
(29). 

We emphasize that the integration over i> and rj effectively amounts to the 
replacement of them by their extremum values. 

The resulting action (30) is quadratic in i?^, and can be conveniently 
rewritten as (here Hi = 112 = f^, mi = m2 = m) 

rT rnp' 1 . 

A= dr\ — + -{aiR'^ + 2a2(Rr)-2ci(Rr)-2c2(rr)+a3r'^ + ay}] (34) 
•^0 ji 2 

where we have used = Rij,-\- [j3 — l/2)r^ to express w in terms of i?, r. We 
have introduced the following notations 

ai = j'^ d^^i + z.), as = j'^ d^{^ + (/? - ^fiy) 

a2 = I' d(3{(3 - i)z., a4 = d(5{^ + ifv) (35) 

/•I A 1 

ci = d[3T]v, C2 = T]v{[5 - -)d(5 
13 



and 

plays the role of the dynamical energy density for the string. Since only 
enters the action then it is convenient to integrate over DR taking into account 
boundary conditions R^{T) — Ri_i{0) = u^T, i.e. 

I DR^ I DR 1^2^ d'^\exp{\^ {R^ - u^Jdr) (36) 

We shall systematically disregard in what follows the preexponential factors 
and shall be interested only in the effective action in the exponent. In this way 
integrating out DR as in (36) we obtain 

I DrDuDrjDfid'^Xexpi-A) (37) 

where 

^ = x / dT{2(Xu) H \ [(aacii - ao)^^ + 2(cia2 - C2ai)(rf) 

2 JO ji ai 

+ {a^ai - ciy + 2a2(Ar) - 2ci(Ar) - A^} (38) 

For generality we preserve the explicitly Lorentz- invariant form of the action, 
with the conditions Rq = l,ro = being implied. It apparently amounts to 
the replacement in eq. (36) of an ordinary integration dXo by the functional one 
DXq(t). The expression (38) for the action A formes the basis of our further 
calculations. 



6 A quantization of the pure straight-line string 

In this section the case of the pure straight-line string will be studied. It 
appeares that in this approximation to the full theory one doesn't need to use 
the conditions (Pr) = 0. We shall show that in this case the reparametrization 
symmetry (B.3) dynamically induces this condition. We shall also obtain in 
this section the spectrum of this particular system. 

In Appendix B it is shown that for the straight line string without kinetic 
terms of quarks at the ends the effective action (38) can be reduced to the 
following one 
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G = J DR^ Dr^ Dr]{T, (3) Du{(3) (39) 

T 1 ^ 

exp[- f dr f d(3—[w'^ + {{(Ti>f + 7y^)r^ - 2r]{wr)]\ 



with being independent on r. 

As well as in the previous section it is convenient to introduce instead of 
z/(/3) a new variable, 

z/(/?) = l/z>(/3) (40) 

which plays the role of an effective energy density of the string. 

Let us first consider an integration over Dr]{T, (3). The functions ?7(r, (3) enter 
the action as an integral over j3 with functions y{l3). Only the extremal values 
of the function i^(/?), which are even under the exchange [3—\ — (/? — |) ( due 
to the symmetry under the permutation of the ends of the string) contributes to 
the action. So we shall integrate only over the class of functions i^(C)(C = 
which are even functions of C,. It is convenient to decompose the functions 
t]{t,(3) into orthogonal polinomials Pn(/5) with the weight z/(/3) 

r]ir,(3) = j:Pnif3)Kir) (41) 

n 

I' d(3u{(3)PMPm{(3) = 5mn (42) 

Taking into account the symmetry of z^(/?), one can easily obtain expressions 
for Po(/?) 

PM = N,, Nl = {j\{p)d(3)-' 

Pi(/3) = (/3-i/2)iVi, Ni = {j\{m-^n?dpr' m 

The action (39) can be written in terms of Ni in the following form 

- ^{rr)ki{T) + \^Nl + 2i{\u)] (44) 
with A being replaced by i\. 
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The function /co(t) enters only in the 4th term in this expression and an 
integration over Dko{r) = UfLi dko{ri) {N oo) gives the factor proportional 
to a product of ^-functions; 

oo 

n ^(Ar(r,)) (45) 

i=i 

Thus only the components of transverse to the direction of A^, which plays 
the role of the total momentum should be taken into account and we have 

I D'^r6{Xr)exp[-S] J D^r exp[-S] (46) 

After the integration (46) the dependence of 5 on A has a simple form 

T(- + 2i{\u)) (47) 

and the integral over A is saturated in the case T — ^ oo we are interested in by 
an extremum value ^ 

= ij^u^ = P^, (48) 
Integration over Dkn with n > 1 leads effectively to the following expression 

ldH0)DUiru)e.p[-l£ dr[l, + (r^ - + Jy^l (49) 

It is important that we have obtained the following constraints 

(ru) ~ (rP) = (50) 

(rp) = (51) 
where is the total momentum of the string and 



P, = (r, - '^r,)^^ (52) 



is the relative momentum of the string. 

The second constraint (51) means that only transverse components of p enter 
the action. The same constraints appear in the canonical quantization of the 
straight-line string [14] . 

Consider the rest system of the meson = (1,0) and transform the expres- 
sions from the Euclidean to the Minkowski space 

dTE — ^ idTM (53) 
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It follows from eq. (49) that the hamiltonian of the problem is 



ff(P>l = ^{^i^ + -VT^" + ]^} (54) 



— * 

where L = {fx p) is the operator of the angular momentum. 

This hamiltonian does not contain the radial part of the kinetic term so that 
the field is not a dynamical one. Noticing that the integral over has the 
form (A. 13), one concludes that in the effective action only the extremum of 
contributes. So that after solving the equation of motion (for a fixed value of 
orbital momentum) 

_i(l±I)iVf + .V^^ = (55) 
we arrive at the final expression for the hamiltonian as follows 



^("■'^^^j^+H/T^'^^^^ ^^^^ 

The function z>'(/3) has the form which gives the minimum of this hamiltonian 
(for details see ref.[13]) 



d(3 



= ^ 8VI(ITT)a 1 

^ ^1 - 4{I3 - 1)2 ^ ^ 

This solution corresponds to the spectrum of the hamiltonian 

Ef = Mf = 27ra^l{l + 1) (58) 

which agrees with the result obtained for the straight-line string in the canonical 
formalism [14]. 

Expression (57) could be physically interpreted if one notices that a param- 
eter 

v{p) = 2{(3 - 1/2) (59) 

plays the role of the velocity of the corresponding elementary piece of the string. 
In this way formula (57) can be rewritten into more familiar one 

= ^ ^\ (60) 
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where pi = ( ^ — ^)^/^ corresponds to the effective mass density of the string 
in the rest frame. 

To conclude this section let us give the physical interpretation of the con- 
straint (50), which shows that (in the rest system) the relative time is unim- 
portant for the dynamics of the pure straight-line string. The given realization 
of the string doesn't include internal interaction between neighbouring points 
of the string, so that there is no physical exchange along the string. Because of 
this fact there is no need for the introduction of relative time dynamics. 

7 The general case of the QCD string with quarks 

This is the central part of our paper. In this section we shall derive the effective 
hamiltonian for the general case of the straight-line QCD-string with quarks and 
shall find the spectrum of the problem. As it has been done in the previous 
section, we integrate over 77(r, /5), expanding it in orthogonal polinomials -Pn(/5) 

— * 

with weight z^(r, /?) After this integration and integration over A in the cm. 
system u = we obtain for the equal mass case (//i(r) = jJL^ij) = I^{{t)) 

S = ]dri^ + ,ir) + lif^+ (61) 

which will form the basis of our further discussions. 

It should be stressed, that the absence of the dependence on the relative time 
ro(T) ( the so-called instantaneous interaction) doesn't mean of course that this 
effective action is induced by the interaction between quarks and string during 
an infinitely small time interval. 

Actually, we have managed to tranform the integration over ro(r) into the 
one over ii{t). And it is the peculiar property of the considered interaction that 
after this transformation one can (under the approximation discussed above) 
arrive at a local three-dimensional dynamics. 

We shall now consider first the case of heavy masses and then arbitrary 
masses with / increasing from / = to / = oo. 

The limit of nonrelativistic potential dynamics. 
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This is the case when m ^ ^/a, so that 

m ~< /J, >;»< u > 
An extremal equation for /i(r) 



2 ^ 
m r 

= 1 + ^ 



/i2(r) 4 

gives after taking into account a nonrelativistic condition 

the simple solution in the leading order 

//(r) = m 

From the extremal condition for i'{t,/3) 



one obtains in the leading order 

v{j3 , r) = cta/^ 
and finally the action in the same approximation is 

m v3 



S = J dr[2m + —r + crv^] 





4 

as one can expect from the very beginning. 

The general case. The transition from the potential dynamics 
to the string dynamics. 



We start with the case / = and arbitrary masses m 

~ (fx f)^ = 
One obtains from (61) with the help of (66) 



^=K[^ + Mr) + ^r^ + .|f 
{ l^ir) ^ 4 
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In the Minkowski space-time the action (70) yields the Hamiltonian 



H = //(r) + 



-9 1 2 

p + m 



+ cr r 



(71) 




and the extremal condition for /i(r) is 



Therefore we arrive at the following Hamiltonian in the case I = 



(72) 



(73) 



where = (pr)^/r^ = 

This expression is widely used in the context of the so-called " relativistic 
quark model" [11]. As was discussed in [2,7] the eigenvalues of (73) differ 
only a little from the approximate version of this Hamiltonian used by one 
of the authors. There eq.(71) has been used with /j, independent of r. As 
a consequence, the procedure was to find first eigenvalues of (71) E{/i) as a 
function of /i, and then to minimize E{ii) with respect to /i, i.e. to find fi = fio 
from ^ = 0, and to calculate E{ii = hq). As can be seen from Table 4 of [2], 
the eigenvalues £n of (73) and En{iJ.o) differ at most by 5% for lowest states, 
while calculations of E{/j,q) are much easier to do, especially for many-quark 
and gluon states. 

Let us discuss how the potential cr | r | obtained in the rest frame is trans- 
formed under Lorentz boosts. One should keep in mind that it is induced by 
the area law of Willson's loop, which is a Lorentz scalar. Therefore it is not 
difficult to verify that this potential is a Lorentz scalar also and in arbitrary 
frame it can be represented as 



where P is the total 4-momentum of the hadron . 

In the case of small values of I, the string contribution to the kinetic part of 
the action (and to the total orbital momentum) 



can be treated as a perturbation of the Hamiltonian (73) (see [2] for a discus- 
sion). If the Hamiltonian (73) (valid only for low I and strictly speaking for 




(74) 
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/ = ) would be used for calculation of the spectrum for arbitrary values of /, 
one could obtain [11] 



= 27ra(2n, + + 5{nr, I)) (75) 

TT 

where X{nr) 4 when / ^ rir and 6{nr^ I) is the small correction for all values 
of rir and /. This formula gives a good approximation of the spectrum for not 
large /. For the Regge trajectory at large I one would get 

= Sal 

However this approximation gives ^ 25% deflection from the correct value, as 
we shall demonstrate below. 

Now we obtain the Hamiltonian of the system for arbitrary I. Separating 
longitudinal and transverse components of r with respect to r , we obtain 

•2 1 • 

= Uff^^ (f X r)^} (76) 
and the kinetic part of the action (61) can be written as 







For the longitudinal and transverse components of the momentum one gets 
respectively 



Vl^^^ = (f (78) 



Pi - = (f + - V2)V(/3, r)dpf^P^ (79) 







The standard derivation of H from the action in terms of these components 
yields in the Minkowski space-time 



(80) 



:f + /(/? - i/2Yvdp) 
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where 

L'^ = {pxr}^ =plf^ (81) 

This is the resulting Hamiltonian for the QCD straight-Hne string with 
quarks. 

Postponing consideration of this general expression to future papers let us 
now concentrate on the transition of the dynamics from the potential case for 
small / with the Hamiltonian (73) to the case of large I, which we call the string 
dynamics. 

In the limit I oo one can expand the potential part of the Hamiltonian 
(80) 

+f^dPr-) (82) 



2\n/2 + J{l}-l/2yvdt3} J V 
around the extremum in | r | 

t] = + l)/(/i/2 + /(/?- \l2fvdi5) ■ J ^dpy/' (83) 
Keeping only quadratic terms in (r — r/) one gets instead of (80) 

,2 



H = y^i^^y^ + Mr) + Ji^d^+ (84) 
V/2 + /(/?- 1/2) W +"^7 ^^^^1^1 



Let us show that in the limit I — ^ oo the dynamical masses /i and i/ satisfy the 
condition 

< II >« > (85) 

In this case expanding eq.(84) in /i/z/, we shall obtain the string dynamics 
regime and the leading Regge trajectory of the Nambu-Goto form 

Ml 27r(7/ (86) 

Now we prove that the alternative regime 

< > > <iy > (87) 

yields larger mass values and therefore it is energetically disfavoured. From eq. 
(84) one can simply estimate that the regime 

< II »< 1/ > (88) 
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is not possible for / oo at all and we are left with a possibility 

< IJ, >^< V > (89) 

which corresponds to the situation when quarks at the ends carry the fraction of 
the total energy (orbital momentum) comparable with the string contribution. 

In this case one can treat the term I{/3 — \/2)'^vd(3 perturbatively. Starting 
with the value ri 

2 _ 1)_2 fd(3^y2 



rf = C-^<7'^yr^ (90) 
and exploiting a numerically small coefficient 

- iiifdfi = ^ 

One returns to the Hamiltonian (73) which leads to the trajectory 

Sal (91) 

with a larger mass for a given / than in eq.(86). This demonstrates that the 
relativistic potential regime is disfavoured, as compared to the "string" one. 

We come back now to the string dynamics regime (85) where quarks at 
the ends carry only small part of the total energy (orbital momentum) of the 
hadron. In the leading approximation neglecting the dynamics of the longitu- 
dinal components one gets 

The extremal value of z/(r, /3) = z//^'' is r-independent, and we recover the case 
of the pure string (56) 



E, --yz.(/?)^i/3+( ) 



2 7 ' V(^-1/2)M/?W 

and hence 

TT 



(93) 



^(o)(^^ ^ ^ 8a(/(/ + l))VY ^,^^ _ _ ^^2)2^-1/2 (94) 
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The corrections to the eq.(92) are considered in Appendix C. It should be 
stressed that while the qualitative dependence of the dynamical mass < ji > 
on I is reasonable 

<^i>r^l'' a>0 l^oo (95) 
<pL> I <v >^ / ^ oo 

the quantitative results are out of the accuracy of the straight-line approxima- 
tion. 



8 Conclusion 

We have applied in this paper the path integral method to quantize the quark- 
antiquark system interacting nonperturbatively. We have argued that the latter 
leads to the appearence of a minimal string between the quarks at large dis- 
tances, R ^ Tg, where Tg ~ 0.2/m is the correlation length in the vacuum. 

As a result our starting Lagrangian consists of kinetic terms for quarks 
(including the relative time term) and the string part. Using the method of 
auxiliary fields one gets an effective action quadratic in coordinates and its 
derivatives. The auxiliary fields are participating in the final action and the 
integration measure for them is found. We have shown that for the spectrum 
the integration amounts to taking the local extremum of the action in the values 
of auxiliary fields. 

One of the important problems dealt with in the paper is the question of the 
relative time. We have shown that this question is resolved for the case of the 
pure straight-line string due to the reparametrization invariance of the action. 

It leads to the constraint (Pr) = 0. which corresponds to the condition 
i?o = in the rest system. We argued that for a string with quarks the same 
constraint should be used in order to make the straight-line approximation 
for the surface selfconsistent. It corresponds to the choice of quark trajectories 
without backward motion in time and defines the valence quark approximation. 
For this class of trajectories it is possible to eliminate the relative time using 
the reparametrization invariance of the string term resulting in a Hamiltonian 
with variable dynamical masses Hiir). This leads to the transformation of the 
initial quadratic kinetic terms to the square - root type ones. 

Thus we have obtained the three dimensional dynamics with dynamical vari- 
ables connected to the dynamical masses of quarks (//^(r)) and to the effective 
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dynamical string mass density v{t^(3). 

The spectrum depends on the relative role of the quark and string degrees 
of freedom (d.o.f.) which leads to the emergence of two different dynamical 
regimes. We have obtained it analytically in two limiting cases, (i) States with 
the total angular momentum / = 0. Here in the kinetic part of the action only 
the quark degrees of freedom contribute while the string provides only the inert 
part, namely, the exactly linear potential. The Hamiltonian appears to be equal 
to the so-called "relativistic quark model" [11] one and eigenvalues coincide 
within the accuracy of our model with those of the proper-time Hamiltonian 
[2,7] where do not depend on r. This potential-like relativistic regime (when 
the string carries only a small part of the total orbital momentum / of the 
hadron) is valid up to moderate values of / smoothly joining the string-like 
regime at large /. (ii) For the states with / ^ 1 the string degrees of freedom 
dominate. The string carries not only energy, but orbital momentum and cannot 
be reduced to the potential term only. The effective Hamiltonian for radial 
degrees of freedom is derived and we show that it yields only small corrections 
to the string levels. 

It is gratifying to note, that the Regge slope of our trajectories is asymptot- 
ically (27rcr)^\ the usual string slope, and it differs by ~ 25% from what one 
would get from the potential regime (i). The latter slope is (8o")^\ the same 
as was obtained long ago in the so-called relativistic quark potential model [11] 
and in [2,7]. It should be stressed, that the string slope is energetically favorable 
as compared to the potential one. 

The case of the pure straight-line string (without quarks at its ends) has 
been treated in detail in Section 6. Results here coincide with those obtained 
by the canonical quantization method [14]. This case corresponds to the limit 
(ii) of / ^ oo of the previous one. 

One should note also, that our results are in an agreement with numerical 
quantization of the same quark-string system, done in ref.[21]. There the in- 
stantaneous dynamics has been assumed from the beginning and quasiclassical 
quantization was performed numerically. The authors [21] also have observed 
two limiting regimes; that of / = and / ^ 1. 
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Appendix A 



A particle Green's function without backward motion 

We shall obtain below the Green's function of a particle in the external 
field under the condition > for its motion. This condition implies 

that the only trajectories taken into account are those without additional pair 
creation, corresponding to the backward in the time motion ^ < 0. 

To illustrate the idea we first consider the case of a free particle subject to 
the same condition 

^ > (A.l) 

The standard form of the free particle Green's function in the Minkowski 
space 

-^Dz^i-f) X (A.2) 

- ^l~s)] = d~s[exp[—{-Y -{x- yf~s)]] 

contains the summation over all trajectories with any sign of In the mo- 
mentum space G has the form 



G{pi , P2) ~ 5'^{pi-P2)^2— — 2 (^-3) 

One can separate in (A.2) the sum over a class of trajectories without back- 
ward motion (A.l) as follows. For each trajectory of this class one makes in a 
unique way the change of integration parameter ^7 by ^ 

di = ^ (A.4) 
and the exponent in (A.2) can be rewritten as 

exv[-^ j ^(^ - im ^ H J^^(^ + «i„(l - f (^„)))] = (A.5) 

^ Z S ^ Z SZq 



T 2 

dT{ ^. + sio(r)(l - z (r)))] 

(t) 
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with 



T = xo - = zo{T) - zo{0) (A.6) 

In the class of trajectories with property (A.l) the transition from the inte- 
gration over dsDzo to (i//(r) with 

//(r) = sio (A.7) 
has a nonsingular Jacobean, well known in the string theory [12] 

Dfi^{r) ~ exp[-i / ^JfI^{r)dr]dsDzo{r) (A.8) 







where ^ ~ A is the ultraviolet cut-off parameter. We note that s plays role of 
a collective coordinate for a set {/i(T)}, since from (A.7) one has 



§=yidTfiiT) (A.9) 



Combining (A. 4- A. 8) one obtains for the Green's function G in the class of 
trajectories (A.l) the following expression 

G = l Dz{T)Df^\T)exp[-i / y(^ + Kr)il - i^(r)))] (A.IO) 







where a proper rescaling of the mass and z{t) is made. Note that in this 
representation the extremum value of //(r) is 



Ij(t) = mVl - ^(r) (A.ll) 

which makes it difficult to treat the case m = and trajectories with z > 1. 

Hence it is convenient to use the canonical, path integral form where one 
gets 

G = J DziT)DpiT)df,iT)exp[i Jip'z- li^^^ + Kr)})dT] (A.12) 







It is important to define the correct integration measure. Having in mind 
that the following equality holds true 
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+ 0O 7, 

/^e:.p[-|(t + l/t)] = 2(^)i/V« (A.13) 
for I arga |< 7r/2, we are to choose (taking into account eq. (A. 8)) 

- n.^. (A.14) 

After the integration over ij.{t) one obtains for G 

T 

G = I Dz{T)Dpexp[i J drip'z - + w?)] (A.15) 



where the exponent is given by the extremum value of //(r): 



//(r) = yjfir) + w? (A.16) 

The expression (A.15) is the usual canonical representation for the quantum 
mechanical Green function with the Hamiltonian 



H = + m? 

In the momentum representation one has instead of eq.(A.3) 



G{pi , P2) ^ S'ipi - n){Ei - Vpf + m2)-i (A.17) 

Thus the contribution of the trajectories without backward motion is equiv- 
alent (for a free particle) to the separation of the positive frequency part from 
the relativistic propagator. 

We take now the case of a particle in the external field A^. In an analogous 
way from the standard form of the Green's function in the external field 

G(x, y\A^)^ J %Dz^exp[-i j d-i{]-C^ - si^) - gz^A^')] (A.18) 







one obtains in the class of trajectories (A.l) the following Green function in 
external field 



T 1 2 

G^nv{x,y I ^^) ~y Dii{j)Dzi,{n)exp[-i j dT{-{j^ + - z ))-gA^,z^ 

(A.19) 

28 



In the same way as it has been done for the free particle to transform 
eq.(A.lO) into eq.(A.15) , one can prove that the expression (A. 19) leads to 
usual relativistic Hamiltonian of the particle in external field 

H{p ,x A^) = -gAo{x, r) + ^J {p + gA{x,r)y + (A.20) 
Appendix B 

The auxiliary field formalism 

To develop a procedure of quantization of (29) we use the auxiliary fields 
formalism, as is usually done in the string theory [3,12]. 
Let us rewrite (29) as 

G = J Dr DRDfj, exp[-K' - K'] Dhabexp[-ao J Vhd^i] ■ (B.l) 
• didaw^dbwi' - habiO) = jDrDR Dhab I D\^^exp[-(j^ j y/hd^^] ■ 

—ioo 



■exp[-i- j VhX'^habd^i]exp[- j VhX'^daW^dbW^d^i]exp[-K - K] 

where 

d^$, = d-fd(5 , 6 = 7, 6 = , h = deth. (B.2) 

In the pure string case when the kinetic terms are absent the action (B.l) is 
invariant under reparametrization 7 — > /(7, /?) 

w;,(7,/3) ^ ^,(/(7,/3),/3) , hn{l,P) ^ (|^)%i(/, /?), (B.3) 
huil. P) ^ {^)hi2{f, /?) , /i22(7, ^ h22{f, (5) 
with the function f{j, (3) satisfying the conditions 

/(0,/3) = ,/(l,/?) = l ,^4^>0 (B.4) 
It appears to be convenient to decompose [12] 

y\o = c^m'^'io ^ r\o (B.5) 
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with 

r%b = ,h^'=ih-Y (B.6) 
The integral (B.l) takes the form 

G = j Dr DR Dhabexp[-K' - K'] J Da{i)Dr\i) (B.7) 

exp[- J{ao-2a{C))Vhd^C]- 

We will show now that in the continuum limit a{^) and f"'^'{^) can be replaced 
by their mean values 

<a(0>^« ,<r'(O>^0 (B.8) 

Equation (B.8) reflects the fact, that a{^) is a scalar, while /"^(O is a trace- 
less tenzor. 

One can write in the form 

a(0=<«(0>(l + KO) (B.9) 
Since a{^) is a scalar, one can expand 

<a(0 >=a + cii?i(0 + ... (B.IO) 

where -Ri(0 is the scalar curvature for the metrics hab', is some constant with 
a magnitude much larger (as we will show below) than a characteristic value of 
i?i(^), so that one can neglect all terms except the first one. 

To simplify the problem we neglect the dependence of hab on ^ and put 
habiO — ^ab- In this casc it is sufficient to consider a model problem, corre- 
sponding to Eq. (B.7) 

J Da{^) DR{-f) Drin) exp{- J d-fa{-f)R^)exp{- J a{-f){f^ + r^)d-f) (B.ll) 



We need the effective action of a{j). The Gaussian integration over DR{j) 
yields 
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The integral over dk is linearly divergent and should be cut off at A; ~ A ^ 
oo. Therefore its singularities in the variable q lie at a distance ~ A 

and the induced action 

Wr^aJ h^{n)d-i (B.14) 

suppresses fluctuations of 6(7). 

The r(7) contribution to the effective action can be obtained with the help 
of transformation 

where f" — f = 0, /(O) = /(I) , and we obtain 



Sr = } dja{j){r' + = / ^(7)i'(7) + / dj{aij)^if'f J ^dj')} (B.16) 
7 

Note that the Jacobian 

does not depend on q;(7). 

As in the previous case one obtains a contribution to the effective action 
from the first term on the r.h.s. of (A. 16) 

Wr^A Jb\'y)d'y (B.18) 

which damps again fluctuations of /(7), and therefore the second term in ex- 
pression (A. 16) is unimportant. 

Thus the field 0(7) is connected in fact with two free fields R{j) and ^(7). 
In this case q; ~ A ^ 00 as shown in [12]. 

In the same way one can show that the effective action for the fields fabiO 
also damps fluctuations of fields 



Wr-^Ajfdj (B.19) 
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Therefore we have justified the equahty (B.8). 
After that we obtain the following expression for G: 



G = j Dr DR D/ii D^i2 Dhab [-K' - K'] ■ (B.20) 

■exp[-{ao - 2a) j Vh(f^] ■ exp[-a j y/hN'^da'w^dbW^cfi] 

with the new action which is quadratic in w;^ and contains the new auxiUary 
fields hab- 

Let us first consider the case of the pure string. The invariance (B.3) makes 
it convenient to introduce the new variables /(O, ^(0? separating out 

the collective mode v{[3) and the field /(7,/?), satisfying conditions (B.4) 



and making a simple rescaling of T112 



n = i^ = {TaV{[3))\ ■'^^;^' f (B.21) 



= 1^ = (^4^)(r/7(7,/?)) (B.22) 
where T enters the boundary condition. Taking into account the fact, that 

Dhii Dh22 Dhi2 = Dh Dhu /i22 ^^22 (B.23) 
and using the well known formula [12] 

Dh'' ~ exp[-^^ [ Vhd^^]Di){f3)Df{-f, /3) (B.24) 

where 1/e ~ A is the ultraviolet cut-off scale, we arrive at the following expres- 
sion after change of the integration over d'y by Tdf{'y, (3) = dr 

G = J DfDh22DrDRDr]Dv[(3)- 
const r 

•exp[-{aQ - 2a H )a / h22{0^{l3)dTdl3]- 

e J 

■exp[- j^dTj'^dft^^({^f-2r,(^r) + ({!>af + rfy)] (B.25) 

where trivial rescaling z,z ^ (^)^^^^' (^)^^^^ together with a proper renor- 
malization of mo, s in K has been done. 
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At first we notice tliat tlie action doesn't depend on /(7, (3) wliicli reflects 
the invariance (B.23). So tliat tlie integral over Df{T,(3) can be factored out 
and it is equal to the volume of the reparametrization group. 

In the standart way [12] we have introduced the physical quantity cr, which 
entered our expression (B.21) 

9 / ^ const. 
cr2 = a{aQ -2a+ —^) (B.26) 

In the general case (B.l) the invariance (B.3) is lost so that we are left with 
a dependence of z> on 7, and have to restrict ourself instead of (B.21), (B.22) 
by a simple rescaling 

n^^^ = (TMT,fi)f (B.27) 

h,^ = ^ = (Tri(t,l})) (B.28) 
After that in the same way we arrive at the expression for G 



G = I Dr DRDji Dh22Dr]Di>{T, (3) (B.29) 

xexp[—K' — K']exp[—{ao — 2a)a J h22{i)^{T, (3)dTd(3] 

xexp[- Jdr Jdf^^ii^f - 2rj{^r) + {(Da)' + Tj'y)] 

which differs from (B.45) by three points; (i) changing factor (ctq — 2a + £211^^ 
by (c"o — 2a) (ii) an explicit dependence of i){T^(3) on r (iii) by the presence of 
kinetic terms. 

Finally after gaussian integration over /i22 > we obtain for the case of the 
pure string 

G = JdR^ Dr^ DV{T)dr]{T,(5)exp[-Astr] (B.30) 

where 




(B.31) 



and in the general case (B.l) 

G = j DR^ Dr^ DiJ, Dv[t, (3)df]{T, (3)exp[-Astr - K' - K'] (B.32) 
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Appendix C 



In this Appendix we derive the effective Hamiltonian for the longitudinal 
excitations, calculate the corrections to the string result (85), and consider the 
behaviour of the dynamical mass ^u- as a function of /. 

It is easy to prove that after taking into account the longitudinal dynamics 
one obtains a small correction to (93) of the order of 

Since fi is weakly growing with /, as we shall see below, one concludes that 
at / ^ oo the field z/(r, /?) separates and is governed by its own purely string 
dynamics, only weakly perturbed by the dynamics of r^(T) and //(r). The latter 
are "living" in the external field iy{/3). 

We now compute the longitudinal contribution to the energy of the system 
to the leading order in (1//). To this end one can use the nonperturbed value 
pj^^ (93) and make an expansion in (83) as follows. 



+ /(;5-l/2)V%^^ V(/5-l/2)Vo)d/?^ 

^ ^2(/(/3-l/2)V0)(i/3)3/2^''/ 4^ (/(/3 - 1/2) V0)(i/3)5/2 

Insertion of v^^\l3) from (93) into (C.l) yields 

{2na{l{l + -2f,+ ^-^{ii)\a{l{l + l))i/2)-i/2 (c.2) 

a/TT 

and finally one gets 

^ = + 2(2.a) V2(/(/ + 1)) + (C.3) 

//^ 12v/2 ^3/2^3/2(^ _ ^^-,2 ^ 

^ 4 ^{an{l + 1))V4 + 23/2(/(/ + l))i/4 J ~ 
where H'^^^ is an effective Hamiltonian for the radial excitations of the hadron 
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After insertion in this expression of the extremal value of /i 

l^ = ('^^f'\^f" (C.4) 
we finally obtain for the effective radial Hamiltonial 

H(r) ^ 1 (34/3( 2 ^1/,^^, ^ ^ {^flH-'l\r - T,f) (C.5) 

where the substitution (/(Z + 1))^/^ — ^ / has been made in the limit I ^ oo. 
Let us consider the case of massless current quarks 

m = (C.6) 

Introducing instead of (r — rg) a new dimensionless variable x 

(r-ro) = 32/¥/iO(^)i/2x (C.7) 

we represent eigenvalues of the Hamiltonian (C.5) in the following way 

< = /-'/'°(^)'/^3^/^«(n.) (C.8) 

where ain^^ is the eigenvalues of the new dimensionless Hamiltonian 

ff"-' = J((-£)'/^ + ^') (C.9) 

In order to obtain the approximate value 0(71^.) we consider eq. (C.3) for the 
restricted class of functions \i independent on r. Such procedure in general 
gives accuracy about 5% for low lying states. In this way one can easily get 

aK) = . 3-3/55(n^ + 1)4/5 (c.lO) 

Substituting this expression into eq. (C.8) we arrive at the final expression 
for the total energy Ei^^r of the hadron and have for the mass squared 

Ml^^ = 27ra{l + constl'^'\nr + 1/2)^/"^) (C.ll) 

which is slightly different in the case / — ^ 00 from the pure string result (85). 
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